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1. Introduction 

Production of top quark pairs in hadronic collision, i.e. at the TeVatron and LHC, is and will 
certainly remain an attractive subject of study. From the multitude of problems, two correspond 
to the current forefront of science: top pair production with two additional partons at the next-to- 
leading order of QCD, and top pair production with up to two partons at the next-to-next-to-leading 
order. Both problems have generated interesting advances in theory. Below, we shall summarize 
the current status. 

2. NLO 

In this proceedings contribution, we present the results obtained with the help of the Helac-NLO 
system. This software consists of the following: Helac-Phegas [jl|], Helac-lLoop [^, based 
on the OPP method and CutTools ||3|], and finally Helac-Dipoles iQ]- 

We consider the process pp^ ttbb+X at NLO ||5|, §] at the LHC, i.e. for = 14 TeV. For the 
top-quark mass, renormalized in the on-shell scheme, we take the numerical value = 172.6 GeV. 
All other QCD partons including b quarks are treated as massless particles. All final-state b quarks 
and gluons with pseudorapidity |t]| < 5 are recombined into jets with separation yA0^ + A)? > 
D = 0.8 in the rapidity-azimuthal-angle plane via the IR-safe ^j- -algorithm. Moreover, we impose 
the following additional cuts on the transverse momenta and the rapidity of two recombined b-jets: 
PTfi > 20 GeV, \y},\ < 2.5. The outgoing (anti)top quarks are neither affected by the jet algorithm 
nor by phase-space cuts. 

We consistently use the CTEQ6 set of parton distribution functions (PDFs), i.e. we take 
CTEQ6L1 PDFs with a 1-loop running a, in LO and CTEQ6M PDFs with a 2-loop running 
in NLO, but the suppressed contribution from b quarks in the initial state has been neglected. The 
number of active flavors is A^^ = 5, and the respective QCD parameters are A^^ = 165 MeV and 
ji^MS _ MeV. In the renormalization of the strong coupling constant, the top-quark loop in the 
gluon self-energy is subtracted at zero momentum. In this scheme the running of is generated 
solely by the contributions of the light-quark and gluon loops. By default, we set the renormaliza- 
tion and factorization scales, /Ir and jXf, to the common value /Xq = m^. 

With these parameters and assumptions we |Q obtain the following results for the cross sec- 
tions at LO and NLO, which are in excellent agreement with |^] 

aLO^(LHC,m, = 176.2 GeV,CTEQ6Ll) = 1489.2 ^'^es.sS ^ > (2-1) 

aJ^lf{LUC,mt = 176.2 GeV,CTEQ6M) = 2636 ^70^(27%) ^ ■ (2-2) 

Besides the agreement for the cross sections, we were also interested in distributions. As far as 
internal consistency checks are concerned, we used the independence of the results on the size 
of the subtraction phase space in the dipole formalism. This led to the plots presented in Fig. |l]. 
The notation and the definition of the phase space cut-off parameter a,„ax is given in the original 
publication. 

The most important result of the study are the distributions presented in Fig. |2[ They show 
the size of the differential K factors, as well as the shape changes in the distributions themselves. 
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Figure 1: Distribution of the invariant mass niy^ of the bottom-anti-bottom pair for pp — > ttbb + X at the 
LHCfor different parts of the real radiation contribution with different choices of Oimax, OCmax — 1 and (Xmax = 
0.01. The red solid line corresponds to the sum of all contributions, the blue dashed line represents the dipole 
subtracted real emission, the cyan dot-dashed line corresponds to the sum of the K and P insertion operators, 
and finally the green dotted line represents the I insertion operator The sum of all the contributions for the 
two different choices of Omax is depicted below. 



It should be noted that based on this complete calculation, the NLO QCD corrections to arbitrary 
distributions can be investigated within realistic selection cuts. On the other hand, this calculation 
served as an important test of the system, as well as a test of the correctness of the preceding results 
from [^]. 

3. NNLO 

The current status of fixed order NNLO calculations for top quark pair production in hadronic 
collisions can be summarized as follows. The most advanced evaluations concern the virtual cor- 
rections. In particular, the latter have been derived in the high energy, fixed angle limit in [^]. 
Subsequently, the complete result for the quark annihilation channel has been given numerically 
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Figure 2: Distribution of the invariant mass m^^ of the bottom-anti-bottom pair (a), distribution in the 
transverse momentum pj^^ of the bottom-anti-bottom pair (b), distribution in the rapidity yj^jj of the bottom- 
anti-bottom pair (c) and distribution in the transverse momentum pj^ of the bottom quark (d) for pp — > 
ttbb -\-X at the LHC at LO ( blue dashed line ) and NLO ( red solid line ). All distributions have been obtained 
with (Xmax = 0.01. 



in whereas leading color and fermionic contributions in the same channel are known from |^]. 
Additionally, the complete divergence structure has been given in | IC]. Thus, the missing contribu- 
tions are the finite part of the virtual corrections in the gluon fusion channel, and the real radiation 
contribution. 

While the full calculation is still under way, there has been important progress in the evalua- 
tion of the behavior at threshold. A first, albeit incorrect, attempt at a complete NNLO threshold 
expansion has been presented in [ 1 1 ]. However, it is only in [ 12] that the known contributions have 
been put together with some missing ingredients due to potential interactions between the quarks 
in order to obtain the correct and complete expansion. 

Assuming that we decompose the cross section as follows 
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where /, j denotes the possible initial states (we concentrate only on gluon fusion, gg, and quark 
annihilation, qq), and with 



1 - 4mj /s 



we have (we have put the scales jx = m) 



+ (- 140.368 In2j8 + 32.106 In j8 + 3.95105) 



(2) _ 3.60774 

+910.2221n^j8 - 1315.531n^j8 +592.292 In^jS +528.5571n/3 +C, 



(2) 

qq ' 



(3.2) 



r{2) 



^^'^f^ + ^ ('496.31n2j8 + 321.1371nj8 - 8.62261 
+46081n'^j3 - 1 894.9 lln^jS -9 12.349 In^jS +2456. 741nj8 +cf' 



gg ' 



(3.3) 



Before we discuss the origin of these formulae, let us note that the cross section expansion in the 
limit, where the emitted gluons are soft (up to two at NNLO), but the final state is not at threshold, 
has been given in [13]. 



The results in Eq. (3.3) require the following ingredients: the two-loop soft anomalous dimen- 
sion at threshold [14, 1^, the matching coefficients at the one-loop level in the two different color 
channels [16] and finally the contribution of the potential interactions between the heavy quarks. 
The latter can be derived from threshold expansion of the top quark pair production cross section 
in e^e^ and 77 collisions [|T7]]. For example the contribution for the gluon fusion is the same as 
that in 77 scattering up to some minor modifications. If we take the formula for the R ratio in 7 
collisions, we have 



GQtNcp 1 



1+Cf 



A(i)+Cf 



A(2) 



(3.4) 



where 



: CfAa + CaAna + TrNlAl + TrNhAh, 



/ 5 1 27 2 25 35 4 ^ 2, /.o^ . 

^- = T2^+(-2 + 8^ JT + T^ +T + T92^ + 

^NA = - ln(2i8)) J + - 1^(2/3)) + ^na; 
A.= (-A + ^h.(2/3))^ + 2x.; 

Ah = 2xH. (3.5) 

From this formula we need to first eliminate the effect of the one-loop matching times the one-loop 
potential (Coulomb) contribution. This is contained in Aa in the term proportional to 1 /j3 . Next 
we retain only the terms with enhancements in j3, i.e. l/j8 or logj3. Finally, we need to take 
into account the possible color channels, since in 7 collisions there is only a singlet contribution, 
whereas gluon fusion requires also the octet. This is done by the replacement Cf — )• Cf — Ca/2. 
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Figure 3: Graphs relevant to the potential contributions in the singlet and octet channels discussed in the 
text. Crosses correspond to the singlet or octet color-projection operators. 

The correctness of the replacement can be checked by considering the color structures depicted in 
Fig. I 



Ref. [12] contains a more extensive discussion of the problem together with a proof that no 



other enhancements are present at threshold at this order of the perturbative expansion. 
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